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In [1-5] there was proposed the method of a factorization of PDE. The method is 
based on reduction of complicated systems to more easy ones (for example, due to 
dimension decrease). This concept is proposed in general case for the arbitrary PDE 
systems, and its concrete investigation is developing for the heat equation case. 

There is considered the category of second order parabolic equations posed on ar- 
bitrary manifolds. In this category, for the given nonlinear heat equation we could 
find morphisms from it to other parabolic equations with the same or a smaller 
number of independent variables. This allows to receive some classes of solutions of 
original equation from the class of all solutions of such a reduced equation. Classifi- 
cation of morphisms (with the selection from every equivalence class of the simplest 
"canonical" representatives) is carried out. There are derived the necessary and suf- 
ficient conditions for canonical morphisms of heat equation to the parabolic equation 
on the other manifold. These conditions are formulated in the differential geometry 
language. 

The comparison with invariant solutions classes, obtained by the Lie group meth- 
ods, is carried out. It is proved that discovered solution classes are richer than 
invariant solution classes, even if we find any (including discontinuous) symmetry 
groups of original equation. 

1 General equation category 

Definition 1. Task is a pair A = [Na, Ea) where Na is a set, is a system of 
equations for graph F C = Ma x Ka of a function u : Ma — ^ Ka- 
Let S {A) be a set of all subsets F C Na satisfying Ea- 

Definition 2. We will say that a (ordered) pair of a tasks A = {Na,Ea), 
B = {Nb, Eb) admits a map Fab ■ Na Nb, if for any F C A^b T e S (B) ^ 

Fx'Bir)eS{A). 

Of course, these definitions are rather informal, but they will be correct when we 
define more exactly the notion " system of equations" and the class of assumed sub- 
sets F C Na- Let's consider the general equation category S, whose objects are tasks 
(with some refinement of the sense of the notion "system of equations"), and mor- 
phisms Mor [A, B) are admitted by the pair {A, B) maps with natural composition 
law. 

For the given task A we could define set Mor {A, A) of all morphisms ^4 in a 
framework of some fixed subcategory A of the general equation category (let's call 
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such morphisms and corresponding tasks B "factorization of A"). The tasks, which 
factorize A, are naturally divided into classes of isomorphic tasks, and morphisms 
Mor {A, ■) are divided into equivalence classes. 

The proposed approach is conceptually close to the developed in [6] approach to 
investigation of dynamical and controlled systems. In this approach as morphisms 
of system A to the system B there are considered smooth maps of the phase space 
of system A to the phase space of system B, which transform solutions (phase 
trajectories) of A to the solutions of S. By contrast, in the approach presented 
here, for the class of all solutions of reduced system B there is a correspondent class 
of such solutions of original system A, which graphs could be projected onto the 
space of dependent and independent variables of B; when we pass to the reduced 
system, the number of dependent variables remains the same, and the number of 
independent variables does not increase. Thus the approach proposed is an analog 
to the sub-object notion (in terminology of [6]) with respect to information about 
original system solutions, though it is closed to the factor-object notion with respect 
to relations between original and reduced systems. 

If G is symmetry group of E^, then natural projection p : N ^ N/G is admitted 
by the pair [A, A/G) in the sense of Definition 2, that is our definition is the gener- 
alization of the reduction by the symmetry group. Instead of this general notion of 
the group analysis we base on a more wide notion " a map admitted by the task" . 
We need not require from the group preserving solution of an interesting class (if 
even such a group should exist) to be continuous and be admissible by original sys- 
tem. So we could obtain more common classes of solutions, than classes of invariant 
solutions of Lie group analysis (though our approach is more laborious owing to 
non-linearity of a system for admissible map). Besides, when we factorize original 
system, defined here factorizing map is a more natural object than the group of 
transformations, operating on space of independent and dependent variables of the 
original task. 

2 Category of parabolic equations 

Let us consider subcategory VS of the general equation category, whose objects are 
second type parabolic equations: 

E :ut = Lu, M = Tx X, K = R, 

where L is differential operator, depending on the time t, defined on the connected 
manifold X, which has the following form in any local coordinates (x*) on X: 

Lu — U"^ {t, X, u) Uij + {t, X, u) UiUj + U" {t, X, u)Ui + q {t, x, u) . 

Here a lower index i denotes partial derivative by x\ form b^^ — is positively 
defined, = Morphisms of VS are all smooth maps admitted by VS task 
pairs. Let us describe this morphisms: 
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Theorem 1. Any morphism of the category VS has the form 

(t, X, u) — > {t' {t) , x {t, x) , U {t, X, u)) . (1) 

Set of isomorphisms of the category VS is the set of all one-to-one maps of kind (1). 

Let us consider full subcategory VS' of the category VS, whose objects are equa- 
tions Ut — Lu, where operator L in local coordinates has the following form: 

Lu — {t, x) {a {t, X, u) Uij + c {t, x, u) UiUj) -\- U {t, x, u)ui + q {t, x, u) , 

and all morphisms are inherited from VS. 

Theorem 2. If set of morphisms Mor-pf {A,B) is nonempty and A e VS', then 
B e VS'. 

3 Category of autonomous parabolic equations 

Let's call the map (1) autonomous, if it has the form 

{t,x,u) ^ {t,x' (x) ,u' {x,u)) . (2) 

Let's call a parabolic equation from the category VS', defined on a Riemann manifold 
X, autonomous, if it has the form: 

Ut — Lu — a (x, u) Aii -|- c {x, u) (Vu)^ -|- ^ {x, u) Vu -\- q{x,u) , ^ (•, u) e T*X 

Theorem 3. Let F : A — > S be a morphism of the category VS, F be an au- 
tonomous map, A be an autonomous equation. Then we could endow with Rie- 
mann metric the manifold, on which B is posed, in such a way, that B becomes an 
autonomous equation. 

Let AVS be the subcategory of VS' , objects of which are autonomous parabolic 
equations, and morphisms are autonomous morphisms of the category VS. 

4 Classification of morphisms of nonlinear heat 
equation 

Let's consider nonlinear heat equation A e AVS, posed on some Riemann manifold 
X: 

Ut = a (u) Au + q{u) . (3) 

(note that any equation Ut = a{u)Au + c(m)(V'u)^ -|- q{u)) is isomorphic to some 
equation (3) in AVS). We will investigate set of morphisms Mor [A,VS) and classes 
of solutions of equation A, corresponding these morphisms. 
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Note, that two morphisms F : A ^ B and F' : A ^ B' are called to be 
equivalent if there exists such isomorphism G : B ^ B' that F' = G o F . From 
the point of view of classes of original task solutions obtained from factorization, 
equivalent morphisms have the same value, that is solution classes are the same 
for these morphisms. So it is interesting to select from any equivalence class of 
the simplest (in some sense) morphism, or such morphism for which the factorized 
equation is the simplest. 

When we classify morphisms for the original equation (2), a form of coefficient 
a {u) is important. We will distinguish such variants: 



a{u) — arbitrary function 



a{u) — ao{u — Uq)' 



a{u) — Qoe 



Xu 



a{u) — ao 



The lower the variant is situated on this scheme, the richer a collection of mor- 
phisms is. Note, that similar relation is observed in the group classification of 
nonlinear heat equation [7]. 

Theorem 4. If a 7^ const then for any morphism of equation (3) into the category 
VS there exists an equivalent in VS autonomous morphism (that is morphism of 
the category AVS). 

Let us give a map p : X ^ X' from the manifold X to the manifold X' and 
a differential operator D on X. We will say that D is projected on X', if such a 
differential operator D' on X' exists that the following diagram is commutative: 

G'^iX') G^{X) 



D' 



D 



C^(X') > G°°{X) 

p* 

Theorem 5. Let a 7^ const. For any morphism of the equation A into the 
category V£ there exists an equivalent in VS autonomous morphism {t,x,u) — > 
(t, y {x) ,v {x, u)) A to S e AVS, for which factorized equation B is Vt = a {v) Lv + 
Q{v), operator L is projection onto Y at map x ^ y{x) of the described below 
operator D (note that this condition is limitation on the projection y{x)), where: 

1) If A is arbitrary (not any of the following special form): D = A, v {x, u) = u. 



2) If A is Ut = UqU {Au + Qqu) + qiu up to shift u 



const: Df = (3^-^ (A {(3f) + qo(3f) for some function (3 
(x) u, Q ^ qiv. 
3) If A is Ut 
for some function /3 : X 



u - Uq, X 0, ao,qo,qi 
X ^ R, v{x,u) 



aoe^" {Au + qo)+qi, A ^ 0, ao, qo, qi = const: Df = e^^ (A/ + A^ + go) 
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V {x,u) = u — P (x), Q = qi. 

We will call such morphisms " canonicaF . In the category VS the canonical 
representative in any class of morphisms is defined uniquely up to diffeomorphism 
of manifold Y, and in the category AVE it is defined uniquely up to conformal 
diffeomorphism of Y. 

Further we restrict ourselves by the investigation of the canonical maps for the 
first variant, that is will look for such maps p from the given Ricmann manifold X 
onto arbitrary Riemann manifolds Y, for which Laplacian on X is projected to some 
operator on Y (note that this canonical maps will be canonical for given X in the 
cases (2), (3) too). 

Note that isomorphic autonomous equations B, factorized given A, arc distin- 
guished only by arbitrary transformations v v'{y, v) and has the same projection 
p : X ^ y{x) up to conformal diffeomorphism of Y. Therefore to find such projection 
p : X — > y for canonical morphism is to find all autonomous morphisms from this 
equivalence class. 

5 Factorizing of heat equation in 

Let VAVS be full subcategory of AVS, whose objects are autonomous parabolic 
equations of divergent shape: 



and morphisms are autonomous morphisms of the category AVS. 

Theorem 6. Let X is a connected region of M.^ with Euchdean metric, F is a 
manifold without boundary, A does not have form (2 or 3) from Theorem 5. Then 
p define canonical morphism of A in VAVS iff p is restriction on X of factorization 
under some (may be discontinuous) group G of isometries. 

6 Factorizing with dimension decrease by 1 

Theorem 7. Let A does not have form (2 or 3) from Theorem 5, and (a) p : X — > F 
is a fibering; (b) X and Y are oriented; (c) X is an open domain in complete Riemann 
space X; (d) dimy = dimX — 1. Then p define canonical morphism to VAVS iff 
the following conditions fulfilled: 

a) p is a superposition of maps pi : X ^ Y' and Pq : Y' ^ Y; 



b) pi : X ^ y is a restriction on X of the projection X — > X / Gi, where Gi is 



c) po : Y' ^ Y is isomeric covering (for the metric on Y', inherited from X); 

d) For the vector field rj generating group Gi, the function '& = {r],r]), defined on 
Y', is projectible on Y. 



Ut — c{x,u) ^ div {k{x,u)Vu) + q{x,u), 




some 1-parameter subgroup of group Isom ( X ) of all isometries of X; 
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7 Factorizing with dimension decrease by 1: com- 
parison with group analysis 

As it was shown in the section 5, when we factorizing heat equation in with 
Euchdean metric, then the class of modeled (3D) solutions of A coincides with a 
class of solutions of A, which are invariant under some (maybe discontinuous) group 
of isometries of M^. 

But this results about coincidence of factorizing maps for the heat equation in 
with Euclidean metric with factormaps by symmetry groups (that is isometries 
groups) are accidental. 

At first, projection pq : Y' ^ Y from previous section is not necessarily generated 
by some group of transformation of Y'. 

At second, let even Y' = Y/Gq, where Go is some discrete group of the isometries 
of Y'. The question is: could group Go be lifted to some group of the isometries of 
X, which preserve projection onto Yl 

Let the group Gi be fixed, which satisfies conditions of Theorem 7. We consider 
differential-geometric connection % on a fibering pi : X ^ Y' with the structural 
group Gi, which horizontal planes are orthogonal to Gi orbits. 

Theorem 8 (necessary condition). If a discrete group Go, which operates 
on Y' and satisfies conditions of the Theorem 7, could be hfted to the subgroup of 
Isom(X), then curvature form dx, projected on Y', would be invariant respectively 
Go. 

Lemma 1. x be decomposed on a sum x = Pi*x' + dh, where x' ^ T*Y', 
/i is a function from X io Ti, H. is fiber of pi (that is either M, or circle M mod H, 
where H — const is integral x on a vertical cycle). 

Theorem 9 (necessary and sufficient condition). A discrete group Gq, 
operating on Y' and satisfying conditions of the Theorem 7, could be lifted to the 
subgroup of Isom {X), iff \/g G Gq the form gx! — x! is: 

- Exact, if the fiber of pi is simply connected; 

- Closed with periods, multiply if the fiber of pi is multiply connected. 
Particularly, \l X — R", and Gi is the rotations group, -q — YlT=i ^i^v'i' "t- > 3, 

or Gi is the screw motions group, rj — 8^ + Yl^i (^i^^fi-i ^ 2, then such groups Gq 
exist, which does not lift on X. 

8 Factorizing with dimension decrease 

Let's equip X with connection generated by planes orthogonal to fibers. 

Theorem 10. Let (a) p : X — > y is a fibering; (b) dimF < dimX. Then p 
define canonical morphism to VAV£ iff the following conditions fulfilled: 

1) The fibers of p are parallel; 

2) The transformation of a fiber over an initial point to a fiber over a final point 
changes volumes proportionally when we translate along any curve on Y] 
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3) The holonomy group saves volume on a fiber. 

Moreover, p define canonical morphism to AVE ifi^ conditions (1-2) fulfilled. 

Example 1 (dimX = 4, dimY = 2). 
Let X — {(x, y, z, w)} with the metric 



9i3 



/ 1 \ 

l + a2 + /52 ^ ^ 

q; 10 

V /? 1/ 



, a 



xe" 



xe 



Y — {{x,y)} with the Euclidean metric, p {x,y, z,w) = {x,y). Then map p and 



equation Vt 



+ Vyy are factorization of the equation 



2/2/ 



2q;m, 



2/2 



-2j3uyy, + (l + Q!^) itz2 + '2.aPuzw + (1 + ^iiytt, + (a/?)^ «z + (a/?)^ 

where a = xe^ and (3 = xe^, by the map p : {x,y,z,w) — > {x,y). (The same 
is true for he equations = a{v)Av on F and it^ = a{u)Au on X for arbitrary 
function a, but for simplicity we will write linear equations in examples.) However 
the only transformations X, under which both the last equation and all it's solutions 
projected by p are invariant, are (x, y, z, w) {x, y, w, z) and identity. Moreover, 
another transformation with such properties does not exist even locally (i.e. it 
couldn't be defined in any small neighborhood on X), even if we replace the request 
"to keep the equation invariant" by the request "to be conformal" . 

Example 2 (dimX = 3, dimF = 2). 
Let X = = {(x, y, z)} with the metric 

(l + z"^ z -z 
z 2-1 
-z -1 1 

Y — {{x,y)} with the Euchdean metric. Let's consider group H of isometries X, 
generated by the screw motion {x, y, z) ^ {x + 1, —y, —z) {H is projectible on y), 

X — X I H,Y — Y I H,p{x,y,z) — {x,y). Y \s homeomorphic to the Mobius band 
without a boundary; X is homeomorphic to the torus without a boundary. 

Then map p and equation e'vt ~ {&^Vx)x + {^^'^y)y-i or Vt — v^x + + on F 
are factorizations of the equation 

Ut = MXO; + Uyy + Mo; + "IZU^^^ + 2Uy^ + + 

on X. However the only transformation X, under which both the last equation 
and all projected by p it's solutions are invariant, is identity map. Moreover, there 
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doesn't exist a non-identity conformal transformation X, under which all projected 
by p solutions of the last equation are invariant. 

Example 3 (dimX = 3, dimF = 1). 
Let X — xM"^ — {(x, y,z) : X e'M. mod l,y^z E R}, equipped with the metric 

gij = \ a 1 , a = -e\ (3 = 2y, 

\P 1/ 

F = S"^ = {a; G M mod 1} equipped with the Euclidean metric, p (x, y, z) — x. Then 
map p and equation vt = v^x on Y are factorizations of the equation 



+ (l + oF) Uyy + (l + P^) Uzz + ^aPuyz " 2aUxy — 2f3uxz + {oi(3)y + (q;/?)^ % 

on X. However the only transformation X, under which both the last equation and 
all projected by p it's solutions are invariant, is identity map. 

Example 4 (dimX = 2, dimF = 1). 
Let X = R^/G with the Euclidean metric, when G is the group generated by the 
sliding symmetry respectively the straight line /. The orthogonal projection of X 
onto the mean circumference (image of the line /) define equation Vt = Vyy on /, 
factorized the equation Ut — u^x + Uyy on X. However the only transformation X, 
under which both the last equation and all projected by p it's solutions are invariant, 
is reflection respectively I. 



9 Factorization without dimension decrease 

If dimX = dimF, then p : X ^ Y projected Laplacian iff it is isometric projection 
up to some conformal transformation Y. 

Example 5. Let's manifold X be a plane without 3 points: A (0, 0), B (1, 0) and 
C (0, 2). Let's consider heat equation on X with metric gij — (x) Sif 

X'^ {x)Ut = UU+U22, (4) 

where A (a;) = p {x, A) p {x, B) p {x,C). p is the distance function (in usual plane 
metric). Let Y = X, and map p : X Y is given by the formula y = — 
^-Y^x^ + ix"^, where x, y are considered as points at a complex plane. 

Because of \yx\ — \x{x — 1) {x — 2i)\ — X{x), heat equation Uf — Uu + U22 on 
Y, equipped by Euclidean metric gij — 5ij, is factorisation of the equation (4) on 

o 

the manifold X, which is obtained by deleting of pre-images of images of zeroes of 

o 

A from X. However, there does not exist non-identical transformation of X, under 
which all projected by p solutions of equation (4) are invariant. Moreover, there 
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does not exist a non-identical transformation of any manifold X\ under which an 
equation (4) is invariant, if X' is obtained by deleting an arbitrary discrete set of 
points from X. 

Example 6. Let's consider an equation on X = R^: 



Let g be the transformation of R^/{0}, that maps x & X to the point, obtained 
from X by inversion under the unit circle with a center in an origin and consequent 
reflection under this center. Equation (5) is invariant with respect to g, but g is 
not defined at origin. However the map p : X = ¥'^ onto the projective plane, 
which past together points x and gx aX x ^ Q/is defined on all X and gives smooth 
projection. Then inducing on Y heat equation is factorization of original equation 
on X. 
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